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Abstract 

\Q ■ We study finite energy topological^/ stable static solutions to a global symmetry break- 

ing model in 3 + 1 dimensions described by an isovector scalar field. The basic features 

^ | of two different types of configurations are studied, corresponding to axially symmetric 

multisolitons with topological charge n, and unstable soliton-antisoliton pairs with zero 
D ■ topological charge. 

oo ■ . 

1 Introduction 

> 

\Q • The familiar solitons in 3 + 1 spacetime dimensions are the monopoles il] of the Yang-Mills- 
Higgs (YMH) model and the Skyrmions j2] of the 0(4) nonlinear sigma model. The first of 
these [T] is a solution of a gauged scalar (Higgs) field model, while the second [2 j is not related 
to a gauge field. But this is not a very strict distinction, since it is also possible to find solitons 
of the SO (3) gauged 0(4) nonlinear sigma model [3]. These are all models that pertain to 
physically rather different contexts but are nonetheless closely related inasfar as they are all 
+f • classical field theories whose static energies are bounded from below by topological charges. On 
this rather technical level therefore, one can ask whether there might be an ungauged Higgs 
analogue of the the (ungauged) Skyrme soliton? We refer to such a symmetry-breaking field 
theory, Goldstone model in 3 + 1 dimensions. 

It is known on the other hand that such models do exist [I] not just in 3 + 1 but in all 
dimensions. These are the gauge decoupled versions of the SO(D) gauged Higgs models |H] 
in D + 1 spacetime dimensions. That such models should support solitons follows from the 
simple fact that certain truncated versions of these have such solutions in closed form [6 i . The 
solitons of these models have the salient feature that the asymptotic behaviour of their solitons 
feature the same properties as gauged Higgs models, and hence afford a simple background for 
the study of Dirac equations [Ij in all dimensions. 

While the existence of solitons to generic such Goldstone models were known for sometime, a 
concrete and detailed construction of these has not been presented in the literature to date. This 
is what we propose to do in this paper, and since 3 + 1 is physically the most relevant dimension, 
we have chosen this for our example. The solitons we construct are in a sense alternatives to the 
usual Skyrmions, though we have not pushed this analogy here. Our strategy here is instead 
to expose the generic properties of such solitons. To this end, we construct the topological 
charge- 1 spherically symmetric soliton, the axially symmetric winding number n multisolitons 
(MS) and examine the possiblity of the existence of bound states. We also construct an (axially 
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symmetric) topological charge-0 soliton-antisoliton (SAS) pair, to highlight the analogy of the 
model studied with the usual YMH model. 

In Section 2 we define the flat space energy density functional of the model and the topo- 
logical charge density presenting its lower bound. In Section 3 we present the charge-1 solitons, 
the charge- n MS, and the charge-0 SAS pairs, in successive Subsections respectively, and in 
Section 4 we summarise our results. 

2 The model and the topological charge 

The symmetry breaking model in 3 spatial dimensions, to which we refer as a Goldstone model, 
is described by a scalar isovector field a , a = 1, 2, 3. There is such a hierarchy of models |lj 
that arise from the gauge decoupled limit of the three dimensional SO (3) gauged Higgs model 
descended from the p-th member of the Yang-Mills (YM) hierarchy on R3 x S Ap ~ 3 . Here we 
have chosen the simplest of these, namely that descended from the 2-nd member of the YM 
hierarchy. Using the notation 

# = 0^ a , <% = dftFdfitp , ^ = 0^0^%^, 

with the brackets [ij...] implying total antisymmetrisation, the static energy density functional 
is 

S {P =2) = A (v 2 - |0t) 4 + Ax ( V 2 - |0 6 | 2 ) 2 K| 2 + A 2 K| 2 , (1) 

which implies a total mass M = l/(47r) J SdV. 

All the dimensionless constants Ao, Ai and A2 must be positive if the topological lower 
bound to be introduced below is to be valid. Moreover, any of these constants can also vanish, 
provided that the absence of the corresponding term in (J2J) does not violate the Derrick scaling 
requirement. That soliton solutions to this model exist is obvious since for particular choices of 
these dimensionless constants, explicit solution jE] are known. Pushing our freedom of choosing 
the numerical values of A , Ai and A 2 further, we can add any other positive definite term to ((H) 
multiplying a new dimension dimensionless coefficient, as long as the scaling properties remain 
satisfied. In 3 spatial dimensions, there is one such possible kinetic term which for which a 
canonical momentum field exists, and that is the sextic term 1 . Thus the most general model 
we can consider is the following extension of (JIJ 

e = \v ( Vj \r\) + r 

t and A being dimensionless constants, k is with dimension of length, and V (77, \<fi a \) a generic 
symmetry breaking potential. 

It is perhaps in order to point out that (J2J) is an ad hoc model, rather than a dimensionally 
descended model like (JTJ). Indeed, a sextic term does appear in the next one to (JTJ 

g( P= 3) = A fog-l^y + Ax (V 2 - I0T) 4 ICI 2 + A 2 (v 2 ~ |0T) 2 K| 2 + A 3 |0^| 2 , (3) 

lr The correspoding sextic term in the Skyrme model was considered in [Hj in some detail. 
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descended from p = 3 YM, with which we are not concerned. We will restrict our attention to 
(J2J), and mainly to the particular case where A = k = 0, which like the Skyrme model captures 
the main qualitative features of the soliton. 

The model (j2j) has certain remarkable similarities to the Skyrme [2] model, and other prop- 
erties that differ fundamentally On the similarity side, there is the obvious shared feature of 
the scaling of the distinct terms in both three dimensional models. Also, for the important 
special case with A = k = 0, the Bogomol'nyi equations are overdetermined like in the Skyrme 
model and there exist no solutions saturating the Bogomol'nyi bound. 

On the contrasting side, the order parameter field <j) a here is a relic of a Higgs field and 
has the same dimensions (L^ 1 ) as a connection, and the finite enrgy conditions require the 
symmetry breaking boundary condition 

lim |0 a | = n. (4) 

r— >oo 

For the (more standard) case of multisolitons centred at the origin, the boundary condition 
there is 

lim |0 a | = . (5) 

r— »0 

For the unit charge spherically symmetric soliton, when the system is described by a single 
function h(r), the conditions (JSJ) and (J3J (see ()10|) below) result in the monopolelike asymptotics 
of our solitons, which are qualitatively different from the instanton like asymptotics of the 
Skyrmions 2 . 

It is straightforward to show that (j2J) is bounded from below by the density 
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whose volume integral is the topolgical charge, which is just the winding number. 

Models like (0), (0) and (jSJ) support global monopoles in the sense that their topological 
charges are the winding numbers of the scalar (Higgs) field on the 2-sphere at infinity. These 
solitons however have finite energy, unlike the usual global monopoles. This is due to two 
reasons. Firsly we have included the quartic kinetic term to satisfy the required scaling, and 
secondly, we have employed a non standard quadratic kinetic term which decays fast enough 
to satify finite energy requirements. These two features are guaranteed by the fact that these 
models are dimensionally descended from higher dimensional Yang-Mills models, the latter 

2 The values of the scalar Higgs function of a unit charge monopole on the boundaries [r = 0, r = oo] are 
[h(0) = 0,h(oo) = 1]. The values of the scalar function w(r) of the spherically symmetric Yang Mills (YM) 
instanton on the boundaries [r = 0, r = oo] are [w(0) = ±1, w(oo) — ±1]. The one dimensional reduced action 
density of the scale invariant YM system in Ap dimensions is proportional to the corresponding 0(2p + 1) 
sigma model reduced action in 2p dimensions, described by the function f(r) = arccosw(r). Hence the values 
of the scalar function cos/(r) of the unit charge sigma model soliton on the boundaries [r = 0,r = oo] are 
[cos/(0) = ±1, cos /(oo) = ±1], like an instanton. 



3 



being endowed with the corresponding topological properties, as explained in and references 
therein. 

While we are exclusively concerned with the classical properties of the model (J2J) here, it 
is nevertheless interesting to comment on its possible quantum aspects. It is believed that 
the Skyrme model is a reliable approximate theory of the Nucleons jSj and our model shares 
many similarities with the latter. The main difference between the models is that while the 
Skyrme field is a constrained field and hence the procedure of quantisation must take account 
of the constraint, the order parameter field in our model is unconstrained. On the other hand, 
the quadratic kinetic term in (J2J) is quite unconventional, featuring the field dependent factor 
(rf — |0 fe | 2 ) . The effect of the latter is to prevent the definition of a propagator, this difficulty of 
quantisation replacing the constraint problem in the Skyrme model. But it is also well known 
that the most efficient practical method of quantisating Skyrme theory is that of collective 
coordinate quantisation, employed in 0. This method applies is equally well to the system (J2J), 
subject to making the essential change of field coordinates 

r a being the three Pauli matrices. 

3 The solitons 

Our aim in this Section is to demonstrate the close similarities of the solitonic solutions in this 
model, with the various monopole (and dipole) solutions of the YMH model. The Section is 
divided into three Subsections. The first involves the charge- 1 soliton analogous to the 't Hooft- 
Polyakov monopole 1J. The second is concerned with axially symmetric solutions satisfying 
standard boundary conditions, with arbitrary winding number n, namely the multisolitons (MS) 
of this model. A question of interest raised in this Subsection is that of the mutual attraction 
or repulsion two 1-solitons. In the last Subsection we impose those boundary conditions on the 
axially symmetric fields, which result in zero charge (unstable) solutions representing soliton- 
antisoliton (SAS) pairs situated on the symmetry axis. These boundary conditions are those 
employed in for the corresponding SO (3) YMH model. 

3.1 Charge- 1 soliton: spherically symmetric 
Subjecting (J2J) to spherical symmetry via the Ansatz 

<j) a = r) h{r) x a , (7) 

and taking into account the factor r 2 in the volume element, the reduced one dimensional 
energy functional, after some rescalings, is 

E = \r]- 2 r 2 V(r},h) +rrj' 
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Figure 1. The profiles of the function h(r) and the mass-energy E(r) of typical spherically symmetric 
solutions are shown for several values of k. 



Mostly, we will retain only the terms multiplying r. In that case the lower bound on the 
integral of (jHJ) is 

Substituting into the limits of the definite integral Q , the asymptotic values 

lirn/t(r) = , lim h(r) = 1 

r—*0 r— >oo 



(10) 



following from (JSJ) and (jlji. one finds 

which was verified numerically. The solutions of the field equation can be constructed numer- 
ically. We follow the usual approach and, by using a standard ordinary differential equation 
solver, we evaluate the initial condition 

at r = 10~ 6 for global tolerance 10 -14 , adjusting for fixed shooting parameter and integrating 
towards r —>■ oo. The behaviour of finite energy solutions as r — ► oo is 

^l + --"-^-^ + 0(l/r 6 ), 

where c is a free parameter. For all considered cases, solutions with the correct asymptotics 
occurs only when the first derivative of the scalar function h(r) evaluated at the origin, h'(0) = b, 
takes on a certain value. For example b ~ 0.443613 for a model without a sextic term (« = 0), 
while b ~ 0.367479 for k — 10 (the symmetry breaking potential is vanishing in both cases). 

The profiles of typical solutions are presented in Figure 1 for several values of the parameter 
k and no symmetry breaking potential. The energy functional, as given by (JHJ) is also exhib- 
ited. No multinode radial solutions were found, although we have no analytical argument for 
their absence. However, our preliminary numerical results indicate that, similar to the case 
of monopoles and sphalerons, the gravitating Goldstone model also presents radial excitations 
with an arbitrary number of nodes of the function h(r). A study of these solutions will be 
presented elsewhere. 

3.2 Charge-n multisoliton: axially symmetric MS 

The axially symmetric Ansatz for the scalar field 

a = (<T,0 3 ) 

is 

(f) a = 7] ip^p, z) n a , 3 = 7]ip 2 (p } z) , (11) 
where p 2 = \x a \ 2 = x\ + x% z = denoting Xi = (x a ,xs), and n a is the unit vector 

n a = (cosn0 , sinn0) (12) 

with azimuthal winding n. 

We denote the two functions (<pi, </? 2 ) = <Pa by labeling <p A with A = 1,2. Subjecting (J2J) to 
(fTTj) . the A = system reduces to 



2\2 



E = 2vrp^r 1 V(l-| V 9 B | 2 ) 



{\d p <p A \ 2 + \d z <p A \ 2 ) + 



2 

P 



6 



2 4 
~T 2 V 



(e AB d p f A d z ip B ) 2 + 



rupi 



(\d pV > A \ 2 +\d zV > A \ 2 ) 
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which in terms of the more useful variables (r, 9) is 
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{r 2 \d r (p A \ 2 + \d e <p A \ 2 ) 



(14) 



and rescaling r as 



t r = | — — | /' — /■ 



removes the coupling constant infront of the quartic term. 

For the multisoliton (MS) solutions at hand, the boundary values of the functions ip A in the 
r ^> 1 region are 

lim (fi(r, 9) = sin9 , lim y?i( r , 9) = cos 9 , (15) 



while at the origin we find 



^1 |r=0 = <f2\r=0 = 0. 



(16) 



Since our imposition of axial symmetry requires also z — > —z reflection symmetry, the actual 
(numerical) integration need be performed only over the range < 9 < ^. The field equations 
have been solved by imposing the following boundary conditions along the axes 



y?i|e=o = d e ^2\e=Q = 0, d e ipi\ e=7r / 2 = <p 2 \ 



0. 



(17) 



We solve numerically the set of two coupled non-linear elliptic partial differential equations 
arising from the variation of the functional (|14j) . subject to the above boundary conditions, 
employing a compactified radial coordinate x = r/(l + r). To obtain axially symmetric so- 
lutions, we start with the n = 1 solution discussed above as initial guess (corresponding to 
(fx = h(r)sm9, tp 2 = h(r) cos 9) and increase the value of n slowly. The iterations converge, 
and repeating the procedure one obtains in this way solutions for arbitrary n. The physical 
values of n are integers. The typical numerical error for the functions is estimated to be lower 
than 10~ 3 . The numerical calculations for n > 1 were performed with the software package 
CADSOL/FIDISOL, based on the Newton- Raphson method [TS]. In Figure 2 we show the local 
mass-energy as given by (|T4*j) of the k = n = 2 MS solution as function of the coordinates 
z = rcos9 and p = rsin^. In Figure 3 the profiles scalar functions <pi and (p 2 of the same 
solution are shown for several angles as a function of the radial coordinate r. 
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E(r,8) 




Figure 2. A three-dimensional plot of the mass-energy E(r,9) of a n = 2, k = axially sym- 
metric MS solution. 

The analysis was carried out in the first place setting the constant k = and for several 
values of n, which captures the main qualitative properties of the MS. The maximum of the 
mass-energy density (JT3jl moves outwards with increasing n. However, for n > 3, the numer- 
ical errors start to increase, and for some n max the numerical iterations fail to converge. The 
problem resides in the behaviour of the scalar function (p 2 , which for large n, tends to develop 
a discontinuity for some value of the radial coordinate. 

Because of the close analogy between our model and the Skyrme model, it is worthwhile 
checking one of the remakable properties of axially symmetric multi-Skyrmions. The property 
in question is that up to vorticity (= baryon number) n = 4 the energy of the multi-Skyrmion 
is smaller than that of n infinitely separated 1-Skyrmions, i.e. that the MS can be regarded as 
a bound state [TH] . 

We have here checked that starting from n — 2, and up to n = 5, the energies of the ra-MSs 
of our model are greater than that of n 1-solitons. Moreover it turns out that this deficit of 
binding energy increases with increasing n, indicating that none of the MSs in this model can 
be regarded as bound states. For example we have found M(n = 2)/(2M(n = 1)) — 1 = 0.229, 
M(n = 3)/(3M(n = 1)) - 1 = 0.381 while M(n = 4)/(4M(n = 1)) - 1 = 0.492 (where 
M(n = 1) = 1.188). 

It is for this reason that we have introduced the sextic 3 term in ()14j) . to check whether 
its presence may reverse this trend and lead to MS bound states? However, we find that for 
< k < 10, MSs with charges up to 3 the deficit of binding energy persists and increases with 
n, confirming that like charged solitons of this model are mutually repulsive inspite of the 

3 The presence of higher order terms in the (covariant) derivatives of the Higgs field is known to result in the 
mutual attraction of like charged (monopoles) solitons |14l 115) . 
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Figure 3. The profiles of the scalar functions tp\ and if2 are shown for the n = 2, k = axially 
symmetric MS solution. 

model having descended from the gauge decoupling of Higgs models [HJ [T3] supporting mutu- 
ally attracting monopoles of like charges. 

3.3 Charge-0 soliton— antisoliton: axially symmetric SAS 

For simplicity we will restrict to winding number 1 SAS solutions, whence we set n = 1 in the 
Ansatz ([Til. 
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Figure 4. A three-dimensional plot of the mass-energy E(r,9) of a n = 1, k = axially sym- 
metric SAS solution. 

The reduced two dimensional energy density functionals (fT3Jl and (|14|) are unchanged for the 
SAS solutions. Also the boundary conditions (fTTj) arising from the requirements of analyticity 
on the z-axis and at the origin, remain unchanged. What are different between the MS and 
SAS solutions are the boundary conditions at r — > oo and at r = 0. 
For the SAS solutions in the region r> 1, instead of (|15j) we require 

lim <f\(r, 6) = sinm9 , lim (p 1 (r, 6) = cos m6 , (18) 

r— >oo r— >oo 

with m > 2 for SAS chains analogous to the monopole antimonopole chains jlbj . But we are 
here only interested in SAS pairs, therefore we restrict to m = 2 in ()18|). 
For the SAS in the region r < 1, instead of ()16|) we require 

y?i|r=o = , d r tp 2 \ r =o = 0. (19) 

The field equations have been solved by using the same methods employed in the MS case. 
However, the SAS solutions exhibit a very different picture. The energy density e = — T\ 
possesses maxima at z = ±d/2 and a saddle point at the origin, and presents the typical form 
exhibited in the literature on MA solutions [TU]. The modulus of the scalar field \<p\ = y^f + 
possesses always two zeros at ±d/2 on the z— symmetry axis. In Figure 4 we plot the mass- 
energy (fT4"j) and the modulus of the scalar field of a typical m = 2 solution as a function of 
the coordinates p, z, for k = (i.e. no sextic term). This solution has a mass M = 2.588 
which is smaller than that of two (n = 1, m = 1) 1-solitons (M=2.92), similar to the sphaleron 
describing a monopole-antimonopole pair ^Uj ■ 



10 



Figure 5. The modulus of the scalar field \<p 
ally symmetric SAS solution. 

4 Summary 



yVi + V 3 ! * s shown for the n = 1, k = axi- 



We have studied the finite energy tope-logically stable static solutions to a (global) symmetry 
breaking model in 3 + 1 dimensions described by an isovector scalar field. Such models can be 
constructed in arbitrary D + l dimensions since they are the gauge decoupled versions of Higgs 
models in all dimensions. We have chosen here D = 3 examples, since this is the dimensionality 
of most physical interest, like the usual Skyrme model, but very different from the latter in 
many essential respects. 

Two classes of solutions are studied: Axially symmetric multisolitons (MS) with topological 
charge n, and unstable soliton-antisoliton (SAS) pairs with zero topological charge, both with 
finite energies. There are two pertinent questions that arise here. In the case of the MS 
solutions, the question is whether solitons of like charge attract or repel, and it was found that 
they always repel, even when the model is augmented with a sextic kinetic term. In the case 
of the SAS pairs, the question is whether they can support a nonzero angular momentum 4 ? 
This task is deferred to some future work, and presumably it will involve a stationary Q-ball 
like features. 

As a scalar theory supporting soliton solutions in 3 space dimensions this model is like 
the Skyrme model. Unlike the latter however this is a symmetry breaking model, as result 
of which the boundary values of the field are akin to that of a monopole rather than that 
of an instanton as is the case for the Skyrme (nonlinear sigma) model. From the viewpoint 
of physical properties, there is no question that it can be regarded as an alternative for the 

4 The correspoding rotating solutions in YMH model have been studied recently in JJj and [T%] . 
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Skyrme model which is known to give a good description of Nucleons j§] at low energies. This 
can be seen from two clear viewpoints: (i) the fact that Skyrmions are capable of forming bound 
states [13] describing exotic states, while the MSs of our model display the opposite property, 
and (ii) because the Skyrmion can be gauged with the (Maxwell) U(l) field [111120] enabling 
the description of the electromagnetic properties of the Nucleons, while the topological lower 
bound on the energy of our MS is invalidated when the scalar field is gauged with U(l). This is 
because Higgs models, from the p = 2 member of which [5] the present model is extracted, can 
be gauged only with SO(D), or, be completely gauge-decoupled as the model considered here. 
While it is true that a 0(D + 1) sigma model in D dimensions can be gauged with all SO(N) 
with N < D ^1 with its energy bounded from below by a gauge invariant topological charge, 
gauging Higgs models with SO(N) with 1 < N < D causes the collapse of the topological 
lower bound on the energy. 

Technically the properties of the model studied here are analogous to those of the usual 
SU(2) Higgs model with symmetry breaking potential supporting monopoles pQ. Like in that 
case, the Bogomol'nyi bound cannot be saturated, and, like-charged solitons repel in the sense 
that a(n axially symmetric) multisoliton of charge n has higher energy than n infinitely sepa- 
rated 1-solitons. (It is possible that solutions with less than axial symmetry may have lower 
mass than the axially symmetric ones studied here, so the possibility exists that such solutions 
may form bound states, however unlikely.) Another point in this analogy is that the sphaleron 
describing soliton-antisoliton pairs here is lighter than two 1-solitons, just as it is in the Higgs 
case [12]. 
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